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Abstract 

We present a new point of view on the quantization of the massive gravitational field, 
namely we use exclusively the quantum framework of the second quantization. The Hilbert 
space of the many-gravitons system is a Fock space ^^(Hgraviton) where the one-particle 
Hilbert space ^graviton Carries the direct sum of two unitary irreducible representations of 
the Poincare group corresponding to two particles of mass m > and spins 2 and 0, respec- 
tively. This Hilbert space is canonically isomorphic to a space of the type Ker{Q)/Im{Q) 
where Q is a gauge charge defined in an extension of the Hilbert space T^graviton generated 
by the gravitational field h^^^y and some ghosts fields , (which are vector Fermi fields) 
and Vfj^ (which are vector field Bose fields.) 

Then we study the self interaction of massive gravity in the causal framework. We 
obtain a solution which goes smoothly to the zero-mass solution of linear quantum gravity 
up to a term depending on the bosonic ghost field. This solution depends on two real 
constants as it should be; these constants are related to the gravitational constant and 
the cosmological constant. In the second order of the perturbation theory we do not need 
a Higgs field, in sharp contrast to Yang-Mills theory. 
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1 Introduction 



The quantization of gravity is an old standing problem of quantum field theory. The solution 
of this problem in full generality is a highly non-trivial problem which seems to be extremely 
complicated. (See however the papers of Ashtekar and collaborators pp, |S21)- In jJEj and 
Plj-P^!, |2n| this problem was addressed for the linear gravitational field of zero mass. Among 
the pioneering works in this approach we mention jSSj, [H], jlE]; QH]; Using the result 

of this analysis many computations have been done in the hterature (see jS], [20] , EE], 

One possible way to perform the quantization of the asymptotic gravitational field is to 
linearize the classical theory of gravitation using the so called Goldberg variables PU] 5 ^H] and 
then to apply straightforward quantization of the resulting free field theory. Because of the 
gauge invariance of the theory (which in this case is the invariance under general coordinates 
transformations) one obtains a constrained system and one tries to use a Bleuler-Gupta type 
formalism, that is to start with a Hilbert space endowed with a sesquilinear non-degenerate 
form and select the physical states as a subspace of the type Qa^ = 0, A = 1, . . . , A^. 

A related idea is to extend the Fock space to an auxiliary Hilbert space Ti^^ including 
some fictitious fields, called ghosts, and construct a gauge charge (i.e. an operator Q verifying 

= 0) such that the physical Hilbert space is TCphys = Ker{Q) / Im{Q) (see for instance [21] 
and references quoted there). As a result of this procedure, it is asserted that the graviton, i.e. 
the elementary quantum particle must be a massless spin 2 particle. The construction of the 
gauge charge relies heavily on classical field theory arguments, because one tries to obtain for 
the quantum gauge transformations expressions of the same type as the general coordinates 
transformations appearing in general relativity. This invariance is then promoted to a quantum 
gauge invariance which should be implemented by the commutator with the gauge charge Q. 

It is an interesting problem to consider the case of massive gravity. This case was analyzed 
many times ago jH], [Ml- In [E] it is argued that even the quantization of the massive spin 2 
field is problematic in the sense that no smooth limit m 1— > exists. Some recent interest on 
this problem exists |2j, |25j, EH and |6j. 

We will show here that one can perform the quantization in such a way that this limit is 
smooth. One finds out that the massive graviton has a scalar partner of the same mass m. The 
construction is done in the spirit of [T^ . 

We also mention that a rigorous construction of the Hilbert space of the many-gravitons 
system is indispensable for the construction of the corresponding S'-matrix in perturbation 
theory in the sense of Bogoliubov. This construction emphasizes the basic role of causality 
in quantum field theory. We obtain a solution for the interaction Lagrangian (the first-order 
chronological product) which goes smoothly for m \ into the solution appearing in [22] • 

The solution we obtain, up to second order of the perturbation theory, coincides with the 
result of the perturbative development of the Einstein-Hilbert Lagrangian with cosmological 
constant, if we make the identification A = 2m^ and use Goldberg variables (see the Conclu- 
sions). We remark that in the second order of the perturbation theory we do not need a Higgs 
field as in the case of Yang- Mills theory. For this reason it seems to be impossible to find our 
massive spin 2 gauge theory by means of the conventional Higgs mechanism. 
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2 The Quantization of the Asymptotic 
Massive Gravitational Field 



One defines the graviton as a certain unitary irreducible representation of tlie Poincare group 
corresponding to zero mass and lielicity 2. In tfie case of massive gravity one sliould use tlie 
representation of positive mass m and spin 2. These representation can be exphcitely described 
using the formahsm of Hilbert space bundles, as presented for instance in (HH , ch. VI.7 thm 6.20. 
Let us denote by HgT^ the one-particle Hilbert space of the graviton of mass m. The ensemble 
of many gravitons is usually described by the associated Fock space J-graviton = JF(+)(Hg'r^) 
where the upper + sign indicates that the gravitons are assumed to be Bosons according the 
the well-known spin-statistics theorem. The Hilbert space J^graviton is not very suitable for the 
construction of the perturbative series of the scattering matrix S in the sense of Bogoliubov. 
The way out is to construct a larger Hilbert space Ti where unphysical degrees of freedom are 
present. In this Hilbert space a (gravitational) gauge charge Q acts which should be chosen 
such that it squares to zero 

= 0; in this case it makes sense to consider the factor space 
Hphys = Ker{Q) / Im{Q) which should be canonically isomorphic to ^graviton- 

Let us describe this construction. We use in this paper the following notations. The upper 
hyperboloid of mass m > is by definition X+ = {p G = m?}; it is a Borel set with 

the Lorentz invariant measure da^{p) = 2^j^- Here the conventions are the following: || -H is the 
Minkowski norm defined by = p-p andp-g is the Minkowski bilinear form: p-q = PoQo—p-^', 
by ?7/ii/ we denote the corresponding flat Minkowski matrix with diagonal elements 1, — 1, — 1, — 1. 
If p G we define r(p) G X+ according to r(p) = {uj{p), p), uj{p) = a/p^ + m^. 

First we consider the zero mass case m = [13], [29] . 

• One generates the Hilbert space H by applying on the vacuum the fields H^^,Up,Up^(^ 
(the rigorous construction is based on the Borchers algebra); these fields are of null mass: 



d'^H^^{x) = d'^Upix) = 



d'^Upix) = ^'^^x) = 



(2.1) 



• H, 



is symmetric and traceless: 



H. 



H% = 



(2.2) 



• The 



field $ is scalar and Hpv,Up,Up have obvious tensor and vector properties 



• The 



fields Hp^y, $ are Bosons and Fermions 



• The 



causal commutation relations of these fields are 



[Hp„{x), Hxu{y)] = - o Vp\Vcjuj + Vpu^VaX - - 





(2.3) 
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and the other (anti) commutators are zero; in particular all Bose fields commute with all 
Fermi fields. Here 

DUx)=Dl:[{x) + Dl^\x) (2.4) 
is the Pauli- Jordan distribution of mass m > and Dm\x) are given by: 

Dt^(^) = ± rf«,t.(p)e^^^-^ (2.5) 

In this Hilbert space there exists a sesqui-linear form (not positively defined) < -, ■ > such 
that we have 

HU = H,u ul = u, ul = -u, <l>t = $ (2.6) 
where by f we mean the adjoint with respect to < •, • > 
The operator Q is well defined through the relations 

gfi = (2.7) 

[Q, H^^] = (^d^u^ + d^u^ - ^rj^^dpU^ [Q, $] = ^^^Up 
{Q, u^} = 0, {Q, u^} = I (^d'H^, + ^9^$^ . 



(2i 



In these conditions one can prove that: 

• The operator Q is well defined; for this one has to check the validity of the Jacoby identity: 

Hx), {/(y), Q}] + {f{y), [Q, b{x)]} = (2.9) 
where b and / are generic Bose (resp. Fermi) fields. 

• The following relations are verified: 

g2 = o (2.10) 

KgQ = QUg, \/g e V. (2.11) 

From ()2.10|) we have 

Im{Q) C Ker{Q) (2.12) 

so it makes sense to consider the factor space Ker{Q)/ Im{Q). One can prove that 
the sesqui-linear form < -, ■ > induces a strictly positively defined scalar product on 
Ker{Q) / Im{Q) and we have a canonical isomorphism Ker{Q) / Im{Q) ~ J^graviUon- 
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The preceding construction presented in detail in justifies the consideration of the auxihary 
Hilbert space ?i as a lay-ground for the perturbation theory. The fields Wp, Wp, $ are called ghost 
fields and the operator Q is the gauge charge. A simplification of the preceding formalism is 
the consideration of the new field 

hf,u = H^, + ^ri^,^ (2.13) 

which is self-adjoint in the sense 

but is not traceless anymore and the causal commutation relations are: 

[hpaix),hxM] = -- iVpxV'^w + VpcoVr^x - VpaVXio) Dq{x - ?/) x 1. (2.15) 

We can easily prove that the preceding definition of the gauge charge is equivalently de- 
scribed by (j2.7j) and: 

[Q, h^J\ = {d^u^ + d^u^ - Tj^^dpUp) 

{Q,m4 = 0, {Q.Up}=id^hp, (2.16) 

so one can consider that the Hilbert space Ti is generated by the fields hpiy,Up,Up with the 
properties described above. 

We now turn to the massive gravitational field. One notices from the very beginning that 
in the case m > the gauge charge defined by ()2.16p does not square to zero anymore. One 
can try to correct this feature as in the case of the massive vector field (see for instance |29J) by 
introducing a new ghost field Vp, which is a vector field. The one modifies the preceding scheme 
as follows: 

• One generates the Hilbert space H by applying on the vacuum the fields hpy,Up,Up,Vp; 
all these fields are of mass m: 

(9^ + m^)hp^{x) = (9^ + m^)up{x) = (9^ + m^)up{x) = {d^ + m^)vp{x) = 

(2.17) 

• hpi, is symmetric: 

hpu = Kp (2.18) 

• The fields hp^,Up,Up,Vp have obvious tensor and vector properties 

• The fields hpy^Vp are Bosons and Up J lip cLFG Fermions 

• The causal commutation relations of these fields are 

[hp„{x), hx^{y)] = -- {ripxri^^ + rjp^ri^x - Vp^Vxu;) Dm{x - y) x 1 

{Up{x),U^{y)} = i Tjp^ Drr,{x - y)l 

[vpix),Vpiy)]='-DUx-y)l (2.19) 

and the other (anti) commutators are zero 
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• In this Hilbert space there exists a sesqui-hnear form (not positively defined) < ■, ■ > such 
that we have 

= V = «M = = (2-20) 

where by f we mean the adjoint with respect to < ■, ■ > 

• The operator Q is well defined through the relations ()2. 71 ) and 

% %TfX 

{Q,u^} = 0, {Q,u^} = i {d'h^. + mv^). (2.21) 

In these conditions one can prove that the operator Q is well defined because of the validity 
of the Jacobi identity ()2.9p and we also have ()2.10|) and ()2.11|) . so again it makes sense to 
consider the factor space Ker{Q) / Im{Q). One can prove in this case also that the sesqui- 
linear form < -, ■ > induces a strictly positively defined scalar product on this factor space. 
However, in this case a modification of the zero mass scheme appears. The one-particle Hilbert 
corresponding to Ker{Q) / Im{Q) is Ht"*'^] © h''^'^! i.e. it describes two particles of mass m, 
one of spins 2 and one of spin 0, respectively. In other words we have Ker{Q)/Im{Q) = 
^gravition ® ^scalar- It scems impossiblc to coustruct a gauge structure such that the scalar 
partner of the graviton is eliminated, so in this paper we will accept that such a particle does 
exists. It remains to be investigated whether the scalar partner of the graviton with a tiny 
mass leads to phenomenological problems. 

Sometimes it is convenient to generalize the expression of the new field fj2.13|) in the sense: 

/^^^ = H,. + \a v,.^ (2.22) 
with a e K*. The causal commutation relations are for this field: 

h^^X^)^ h^liv) = {vpWau + Vpu^V^x - VpaVxu^ Dmix - y) x 1. (2.23) 

We can prove that the definition of the gauge charge is equivalently described by ()2.7j) and: 
[Q. h^^J] = -l + d^u, - Vpudy^ , (2.24) 

{Q, u,} = 0, {g, u,} = I (^d^h^^^ + ^ d,h^-^ + mv)j (2.25) 



and 



here 



[Q.v,] = -'^u,- (2.26) 



f^{a) ^ ^pu j^{a) _ ^2.27) 
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The choice (I2.13|) correspond to a = 1. Let us consider the choice a 
preceding relations for 



become: 



h = h^-^^ 



[Q,h 



u'^p) 1 



and 



\n 1 



-1. Then the 
(2.28) 

(2.29) 
(2.30) 
(2.31) 

(2.32) 



This choice seems to appear naturally in the classical framework of gravity with an non-zero 
cosmo logical constant, if one expands the metric g^i, around Minkowski background in the form 



9tMu = rii^y + 



(see the Conclusions). However, from the quantum point of view the value of a is irrelevant: 
all choices are good for the description of the physical Hilbert space. 

We remark also that the massless limit problem mentioned in [S] has a very simple expla- 
nation according to the preceding observation: in j3] one uses different values of the parameter 
a for the case m = and m > respectively. The correct procedure is to use the same value 
of a in both cases. 

The construction of observables can be done in the usual way. We denote by W the linear 
space of all Wick monomials on the Fock space 7i^^ i.e. containing the fields h^^{x), u^{x), 
u^{x) and f^(a;). If M is such a Wick monomial, we define by gh±{M) the degree in (resp. 
in u^). The total degree of M is 



deg{M) = gh+{M) + gh_{M). 
The ghost number is, by definition, the expression: 

gh{M)=gh^{M)-gh_{M). 
If M G W let us define the operator: 

cIqM =. QM : -(-1)9^W : MQ : 
on monomials M and extend it by linearity to the whole W. Then dqM G W and 

ghi^dgM) = gh{M) - 1. 



(2.33) 



(2.34) 



(2.35) 



(2.36) 



The operator dg : W ^ W is called the gauge variation; the properties of this object are 
summarized in the following relations: 

dl = (2.37) 
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2 

dqUf, = dqUf, = i {d'^h^y + mv^) 

tin 

dqVf, = — —Uf,{x). (2.38) 

dqiMN) = {dQM)N + {-ly^'^^''^ M {dqN) , ^M, N eW. (2.39) 
If O : 7Y^^ — s> H^'^ verifies the condition 

dqO = (2.40) 



then it induces a well defined operator [O] on the factor space Ker{Q)/ Im(Q). 

Moreover, in this case the following formula is true for the matrix elements of the factorized 
operator [O]: 

m,[o]m = {^,o<^). (2.41) 

This kind of observables on the physical space will also be called gauge invariant observables. 
An operator O : Ti.^^ —>■ induces a gauge invariant observables if and only if it verifies: 

dQO\Keri,Q) = 0- (2-42) 

Not all operators verifying the condition ()2.4U|) are interesting. In fact, the operators of the 
type dqO are inducing a null operator on the factor space; explicitly we have: 

[dQO] = 0. (2.43) 

In the framework of perturbative quantum field theory the axiom of factorization in the 
adiabatic limit is: XX 

lirndq dxTr,{xi,...,Xn)g{ex)\j^,n) = 0, Vn G N*. (2.44) 

If infrared divergences cannot be avoided, the one can consider the preceding relation at the 
heuristic level and impose the postulate: 

dqTn{x,,...,Xn)=tJ2gniT!:/i{xi,...,Xn), G N* (2.45) 

1=1 

as it is done in [2H]- In particular we have for n = 1 

dqT,{x)=i^T^{x) (2.46) 

for some Wick polynomials T^{x). The derivation of the most general expression of Ti can be 
done in the original variables H^^, $, -u^, m^, or with /i^"^ for any values of a. If we change the 
fields, we must correspondingly change the expression of the operator dq and the final result 
should be the same. In formulae: 

dqT,{H,,,...)=z^T^^iH,,,...) ^ dqT,{h^^J ,...) = i^T^{h^^;} .) (2.47) 

for any a. 
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3 First Order Gauge Invariance 

In view of the discussion from the preceding Section it is natural to discard from the interaction 
Lagrangian (the first-order chronological product Ti) expressions of the type 

dQB + d^B^, gh{B) = -l,gh{B^) = 0- (3.1) 

we call such an expression a trivial coupling. If the difference of two couplings is a trivial one 
then we call them equivalent. In this Section we prove the following 

Theorem 3.1 Let us consider the most general Wick polynomial T tri-linear in the fields 
H^y, ^,u^,u^,v^ verifying the following conditions: 

UgT = TUg, MgeV 
gh{T) = 

3 < deg{T) < 5 (3.2) 
and the gauge invariance condition \2.4(^ - Then T is equivalent to the following expression: 

T = a T^"'^ + b T^''^ a,beR (3.3) 

where 

-4m{dpV,)u>'u'' + (^^H'^" HppHJ> + ^H^^'H^,^ - - ^u^u^ + j (3.4) 

and 

+m[H^''{dpH^,)vP - 2m''{d,HjvP + H^''{d^^)v, + (9X)wX - l{d^.v,)u^u' 

(^-^H^^-^H^pHf - H^^'upu, + '^-H^^'v^v^j . (3.5) 

The preceding expression has a smooth limit m 0. 

Proof: We make a list of all Wick monomials verifying the conditions (|3.2p . The condition 
of Lorentz invariance depends essentially on the dimension of the space-time. In other dimen- 
sions than 4 the list below changes drastically. Also the fact that H^i, is traceless is useful in 
eliminating many terms. First we have terms which do also appear for the massless case i.e. 
without the ghost field namely: 
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rp{Z) ^ ^(3)^3 

2^(8) _ c^^)^Hf^{df,d''H^^) 

2^(10) _ c(^°)(f^''$)(5"'^)^A.. 
2^(16) ^^(15) 

then we have the terms containing at least one factor namely: 

C/(^) = d^^^Hi'''{d^^)v, 



+d'. 



(12), 



^(13) _ (^(13)$^a(5^5^^/3 

lj{i5) ^ d^^^^v'^v^idadpdX)- 



(3.7) 



We have discarded a lot of terms because up to a total derivatives they are of the type 

already considered. As a general strategy, wc have eliminated all terms with derivatives on 



the ghost fields u^. It is somewhat more complicated to prove that one can make c 



one subtracts a total divergence and redefines cf^ , cf^ , c^^ , and Cg 
'ergence and r 
Also because 



(6) 



(6) 
3 



if 



if one subtract a total 



divergence and redefines cf^ , cf^ . 



(a2 + m5)/, = 0, J -1,2,3 =^ 
{d'h){d,f2)h = \{ml + ml- ml) hhh + \d, [(^VOMa + h{d^f2)h - fMd'^h)] (3.8) 

we can eliminate many terms by subtracting a total divergence and redefining other terms of 
lower canonical dimension. 

Now we can put to zero some of the constants above if we subtract from T a coboundary 
i.e. an expression of the form dqB where we take S to be a Wick polynomial with the following 
properties: 



gh{T) = -1 
2 < deg{T) < 4. 



(3.9) 



We have the following admissible expressions: 

First we have terms which do also appear for the massless case i.e. without the ghost field 
V/^ namely: 



10 



^(3) ^ fe(3)$//M-(9^^^^) 

= 6f )tiX(5p^^'') + + bf^^^rid^u^) + hf e^,p,u>'u''{dPu''); (3.10) 
then we have the terms containing at least one factor namely: 

yd) = fWH^-v^u, 

y(3) ^ f[^)y-(d^y^)u^ + fPv^{d^Va)u^ + V,(d''v^)u^ + fi'^e,,apv''(d''v^)u^. (3.11) 

We have discarded some terms because up to a total derivatives they are of the type already 
considered. Now one can prove that: 

• One can use B^^'> to make C2 \cf \cf\cf^ equal to zero. In this way one redefines 

T(7)^ 7^(8) ^2^(11) ^2^(12)^ 2^(13) _ 

• One can use B^'^^ to make equal to zero. In this way one redefines T^^\T^^^\T^^'^\ 

2^(14) ^2^(15)^ 2^(16) _ 

• One can use B^^^ to make c^^^ equal to zero. In this way one redefines T^^^ — T^^^\ 

2^(14)^ 2^(15) _ 

• One can use S^^^ to make c^^^^ equal to zero. In this way one redefines T^^\t^^^\ 

• One can use B^^^ to make ci'^\ c^'^\ c^l'^\ c^^'^^ equal to zero. In this way one redefines 

• One can use V^^") to make d'^^ equal to zero. In this way one redefines T^^\ U^^\ U^'^\ 

• One can use V^^^ to make o?^^-* equal to zero. In this way one redefines T^^\U^^\U^'^\ 

• One can use to make XX d^i\d'i\ d^ , d^P equal to zero. In this way one redefines 

C/(7) _[/(10)^t;(12)^^(13)_ 

One can count that we are left with 37 coefficients of type c^''^ and 23 coefficients of type 
d^\ So we have 60 free parameters which should be fixed by the condition of gauge invariance. 
Now one considers the coupling 

16 13 

T = ^ T(^) + ^ (3.12) 
i=i j=i 
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and computes the expression dqT. It is convenient to follow the same strategy as above and 
eliminate, up to a divergence, the derivatives on the ghost fields u^. One gets by a straightfor- 
ward but tedious computation the following result: 



dqT = i d^X^" + iu^'Yi, + i Z^ + i Ze 



(3.13) 



where the expressions X'^, do not contain ghost fields and the expressions Zj are tri-linear 
in the ghost fields of canonical dimension j = 4, 6. The explicit expression of X'^ and Zj are not 
important for the moment. The expression of is rather long: we group the various sectors: 



22 



where 



= 
— 



3c« + 771^4') --mdg 



T III- — —TTld^ ^ 

3 (1) ^ ^(2) ^ ^2^(6) ^ 1 ^(3) 



y2 = 



Y^ = 



Y' = 



2 2 2 2 

3c(3) + + 1^2^(9) + ^2^(10) _ 1 2^(14) ^ 1 ^(6)1 ^(g^^) 



(3.14) 



(3.15) 

(3.16) 
(3.17) 
(3.18) 



^ lcf\d,d^H,^WHn + [-2cf^ + - ci 



(a.a''i/^-)(a,i/^,) 



.(6) 



+ C2 



,(7) 



,(6) , J6) (11) 
1 "T Cg C4 



2cf ) + 4^^) 



J6) 



+ -C(«) + C5 



,(11) 



<-6 ^ '"3 



(a^a^i/p,)(a,i/^-) 



+ 



24^) + 4") 



2c^') - c^'^ 



+ 



J6) 



— C 



.(11) 



-ci^'^e^^,{d'^d^H,^){d''H^^) (3.19) 
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2cf ) + cf ) 



pa) 



+ 



-.(6) ^(13) 



pa) 



(3.20) 



.(7) 



^2 ~'~ 2*""' 



(11) 



1 m 1 



+ 



2 ^ 2 



(10) _ 1 (11) _ 1 (11) _ _(16) 



+ 



.(7) 



-c, 



(11) 



-c 



(11) 



e^.p/3(9'^9,$)(9^//^^) (3.21) 



(a,$)(a''a'^i/^,) + 



(a,$)(a,a^i/'"^) 



+ 



-c. 



(7) _ ^(9) _ ,(14) 



y9 



y; 
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((?^$)(a,api/'^0 + 



C3 -l- C4 



yii = :^ 
M 2 



\^12 = _ 
— 



-.(10) _ ,(14) _ ,(16) 



e,.,p{d''^){d,d^Hn (3.22) 

(3.23) 



y 



13 



1 



1 
2' 



,2^(8) 



y 



14 



(10) 



mc 



mc, 



f ) _ irf(4) + 47) _ I^^^(IO)- 

-cf) + - Irndr^' 

ai) _ 3^(3) _ 1^(4) _ 1 410)" 



{d,H,,){d''v'') 
{d,H^){dfv'^) 



.(11) 



— -md 



(10) 



(3.24) 
(3.25) 
(3.26) 

(3.27) 
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+ 
.(11) 



.(11) 



(3) 



-nid[ 



(10) 



(10) 



1 r 
2 



3df ^ + mdf ) 



^15 _ 



.(13) 



(10) 



mc^^^ H — dg"^'' — md\ 



j(3) 



(10) 



+ - 



143)^^(4)^/10)^^^(10) 



(13) 

mc\ 



-mdf^ — md^^^^ 



^16 _ 
M — 



mc. 



+m 



.(14) 



- -d 



1 r 
+ 2 

(7) J 



d 



(12) 



yi7^ _ 
M — 



mc 



(16) _ 1^(4) + 1^(8) + 1^(8) 1^^(13) 

4 2 ^ 2 ^ 2 



Y. 
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^(1) _ mdP - mdf + -mci(9) - Im^d^^ 



+ 



rf(l) _ ^47) _ ^^(8) _ 1 2^J10) 



1^(1) + ^(2) + ^47) ^ ^^8) ^ 1 ^(9) ^ ^2^(10) 



— m 



4^^ + 4"'' — '"^4 

2 



/(8) 
<'4 



(10) 



Y. 
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14'^ + 14^) - 4'^ 



-/i' + d? 



df + df^ 



14 



Y, 
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(12) 



(3.34) 



{d^v.WdPvp) 



+ \dT\d,v,Wd^Vp) + ]^df\d''v''){d,d,v,) 



+ 



1 



y22 ^ _id(i3)^«(a^a„a,.;^). 



(3.35) 



(3.36) 



We now impose the gauge invariance condition ()2.46p . It is sufficient to take a tri- linear 
expression in the fields i/^^jy, $, m^^, m^^, verifying the following conditions: 



If we define 



then we get from (|2.46p and (|3.13p : 



gh{T) = 1 
3 < degiT^") < 5. 



(3.37) 
(3.38) 

(3.39) 



The generic form for is 

T'^ = T^'" + (dpu,) T^^'P + {dpd^u,) T'^'^P'^ + dou^u^'u, + die^'P'^u.UpU^ + (3.40) 

where the expressions yM^"^ ^mj^p^ t^m'^pc g^^^^ bi-linear in the fields H^i,,^, A^, the expression S'^ 

is tri-linear in the ghost fields and of canonical dimension 5 and rfo, di are constants. We have 
by direct computation: 

+ {d^dpd^u,) T^'^P'^ + do [{dpu^)u'u, + u^idpu^u, + uV{dpU,)] 

+rfie'^^'"^ [2{dpU,)upU^ + u^UpidpU^)] + dpS^'. (3.41) 



Let us write 



rpfll/pcr rp/MUpa _|_ rppupa 



(3.42) 
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where the tensor T^^^" does not contain terms with the factors -qi^P ^ rj'^'^ , t]'^'^ and collects 
all terms containing at least one of these factors. The terms with the factor rj'^'' can be discarded 
if we redefine T^*^. It follows that the generic form is 

T^^- = 1 (r^^^rr + v'^^T^') ■ (3.43) 

We also write 

jnpup _ rppup _|_ ^pprpu ^2 

where T^'^'' collects all terms without the factor 7]^^. 
Then we obtain from ()3.4ip 

a^T'^ = {d^T^"" - m^T") + {dpu,) {TP" + d^T^^^ + dfT" - m^TP") 

+d^eP'""' [2{d^u,)upU, + u,Up{d^u^)] + dpSP. (3.45) 
It follows that the equation ()3.39|1 is equivalent to the following system: 

TP'' = -dpTl"'P - dPT" + m^TP" 

{d,SP), = Zfi. (3.46) 

Here we have written 

where the two pieces have the following symmetry properties: 

j.pup^^j,pvp Ve = ±. (3.48) 

The expression {d^SP')j contains the terms of canonical dimension j = 4, 6 of df^S^] the 
expression ((9^S'^)4 is proportional to the mass, so is zero in the massless case. 
From the first four equations of the system ()3.46|) we obtain 

Y'' = {d^ + m')(9^Tr - T'^), (3.49) 

so Y^ has the generic form 

Y^ = ]^{d^ + m')Zp (3.50) 
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where is a Wick polynomial bilinear in the fields H^j^^, $, and verifies the following con- 
ditions: 

gh{T) = 

2 < deg{T^') < 3. (3.51) 

The generic form of is 

= fMd^^) + f2H^''{d^H^p) + hHd^H^,) + hH,,{d''<l>) + hH,,{dpH^P) 

+fiiv''{d^v^) + fuv'^idavn + /i3 W^^l^'^'')- (3-52) 
The basic equation (j3.5(J|) becomes now 

If we substitute here the expressions fj3.14|) - ()3.36|) we get the a system of equations of 
78 equations for 73 the unknowns c^^\ fj (37 + 23 + 13 = 73). One can solve this system 
explicitly, the only non-zero coefficients are: 

c^"^-* = -rn^a m?h c^^-* = —m?a c^^'' = — -m^a c*-^^ = —m?h c*-^-* = —ni^a. 

3 3 2 4 

cf = -2a cf = -Aa = 2a c^^) = 4a c^^) = a 
c(^°) = a 4^^)= 4a cf) = 4a cf ^ = 4a 4^^^ = 2a 

= -m% d(^) = m^a = mb = -2m6 rf^^^ = mb 
d^'''^ = mb = —Ama — -mb d^p = mb d^^^ = mb d^2^ = mb df^ = —-mb 



= 2mb dS'°^ = 4a - 6 4'°^ = 2& 4'°^ = -2& (3.54) 



and 



/i = -^a /2 = -a /4 = -2a /s = 4a /g = 4a; 

/g = -mb fio = b /n = 2a - ^6 fu = -b. (3.55) 

here a, 6 G M are arbitrary parameters. Now we can determine the expression Zq i.e. the terms 
tri-linear in the ghost fields and of canonical dimension 5: 

Ze = 4a u"(a^M^)(9^9^M„) (3.56) 
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and the last equation of the system (j3.46j) admits a (non-unique) solution. One can choose the 
"minimal" solution: 

= 2a [u^{d^'Up){d''uP) + u^{d^'Up){dPr) - u^Up{d''d''uP)] . (3.57) 
Finally we consider the fifth equation (j3.46|) . First we have the explicit expression: 



Z^ = — 
2 

1 

+ 2 



J- 2 (11) 

-m Cn 

4 



-L 2 (13) 

-m Co 
2 ^ 



(7) 



+- 



(3.58) 



It follows easily that we can take rfi = in the expression ()3.40|) of T''. Then we get from 
the fifth equation (|3.46|) : 

3 

d^ = —2m^a — -m^b, (3.59) 



so we get the result from the statement. ■ 

We remark that the basic idea of solving the gauge invariance problem was to use the 
equation ()3.50|) instead of the original gauge invariance condition 1)2.46^ : the equation ()3.50|) is 
simpler because we do not need an ansatz for T'^; only an ansatz for Z'^ is necessary. 

The solution T^°-^ gives in the massless limit the usual gravity theory [221 plus the new term 
4:H^^'^{d^Vp){d^vP). The usual choice is a = — |. 

One can re-express using the variables from ()2.13p : we have 



Proposition 3.2 In the variables 

h^u = + ^r/^^$ h = W^ = 2$ 
the expression T^"^ from the preceding theorem is equivalent to: 

-2/1'^- (9, V)!^"^) + '^h^^''{d.hpp){d'^h-^P) 
+Aid^hnuPidpU,) - Ah'^'idX)id.Up) + Ah^''idpVp)id,vP) 

4 „, . 1 



(3.60) 



In these conditions one can take in 

T^-) = u''[-id,h){d,h) + 2{dphp,)id,hn - Aidphnid.h^.) 



(3.61) 



-{d.hWh) - {d^hp^Whn - 2{dpK^){d^h''P) + 2{d,Vp){d''vP) 
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+m'u^ (h'^^Kp - ]^h^\ . (3.62) 



Proof: It is convenient to start from the expression T above and make the substitution 
H^^, = h^jj — h = 2$. Then one makes the transformations described at the begin- 

ning of the proof of the preceding theorem, namely we get rid of derivatives appearing on Up 
subtraction total divergences. We also note that the fourth and the fifth terms from the expres- 
sion of T can be eliminated if we use the identity ()3.8|) but have been included such that the 
linear approximation of the Hilbert Lagrangian is reproduced [21]. We obtain the expression 
from the statement of the theorem. 

The computation of Tjf"^ is not very difficult and provides an impressive check that the 
preceding computations are right. It is not hard to obtain by direct computation that 

dgT = tlu^t'' + + id,dpUp)t''''P + s] (3.63) 

where the expressions t^^t^^^t^^P are bi-linear in the fields h^^^v^ and s is tri-linear in the 
ghost fields. We note that we have a certain freedom in choosing the expression t^^P if we 
do not impose symmetry in z/, p. Now, as in the proof of the preceding theorem, one makes 
"integrations by parts" and rewrites the preceding expressions as follows: 

dqT = liu^'y^ + d'^x^ + s) (3.64) 

where 

yp = tp- + d^dHp^p; (3.65) 

and 



Xp = u^Xpi, + {dPu") 
•^pu — tufi d^ti^pp Xpvp = tupp- (3.66) 

By direct computation we can prove that 

Vp = (3.67) 

It remains to prove that the expression from the statement is such that s = d^'Sp, and the 
expression dqT above is exhibited as the total divergence from the statement. ■ 

In the same way one can also re-express T*^^^ using the variables from (3.60); we have 

Proposition 3.3 The expression T^^^ from the preceding theorem is equivalent to: 

T = -h'^-{d,Vp)id,vP) + 2h>^''idpv,){d''vp) - ^v{dpvPf 
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+m 



In these conditions one can take in \2.4(^ 

T^^^ = u^{{d,v,){d,v'') - {d,v,){dX) - \{d,v.){dpv'') 



+m 



-{dphpy)vP + {d^hpp)v'' + ]^{dpKp)vP - ^{dp,h)v^ 
-]^{duh)vp + hi{dpv^) - ]^Kp{dpVp) 
(^-^VpV^ + hpph'P - ^hp^h^ | 
+u,{-l{d.VpWvn + lidpVpf - lid.vp){d^vn 
-h'-idpv^) + \h{dpvf) + \{dph)vP - ^{dPhp^X 
W(-^V^"' + + ^v-^vj^^+mid^unhppv, 
-lidpu,) [v'id^vp) - Vpid^vn -m{h''Pvp - i/it;'^)] —m'upU.ur 



+m 



(3.68) 



(3.69) 



20 



4 Second Order Gauge Invariance 



In this Section we consider the second order gauge invariance. For this we must construct 
the chronological products T{x,y) and T^{x,y) such that the identity ()3.39|) is verified. The 
construction procedure is well-known: one computes first the corresponding causal commuta- 
tors \T{x),T{y)] and [T^(x), T(?/)] and makes the substitution Dm{x — y) ^ D^{x — y) i.e. 
one substitutes the causal Pauli- Jordan distribution by the corresponding Feynman propagator 
and similar substitutions for the loop graphs; one obtains the expressions T^{x, y) and T^{x, y) 
which verify all Bogoliubov axioms but might spoil second order gauge invariance. To restore 
it we must annihilate some anomalies and make finite renormalizations. These finite renormal- 
izations must also preserve the power counting theorem which in this case says [211 that the 
expressions T{x,y) and T^{x,y) should be of the form 

T(x,y) = ^t,(x-y)iy,(x,y) (4.1) 

j 

where Wj are Wick polynomials and tj are distributions such that one has 

Lu{tj) + deg{Wj) < 6. (4.2) 
The origin of the anomalies is explained in ^2]; [21] ■ One starts from the identity 



dQ[nx),ny)] = ^^[T''{x),T{y)]+^^^[T{x),T>^{y)] (4.3) 

which follows from first order gauge invariance. If one substitutes expressions of the type ()4.H) 
for the causal commutators 

Dix,y) = J2djix-y)W,ix,y) (4.4) 
j 

then the preceding identity reduces to some identities verified by the distributions dj. When we 
make the causal splitting of the these distributions, preserving the degree of singularity, some 
of these identities are lost and we get anomalies. From tree Feynman graphs we get the identity 

{d^ + m^)D^ = (}, (4.5) 

which cannot be split causally preserving the degree of singularity; indeed it is well known that 

{d^ + m^)D^ = 6{x-y). (4.6) 

From loop graphs we get the identities 

d,D^\,n.2 = -<D:^2,m, + D^^^^^^ (4.8) 
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d,D^Z,m, = -^2^:;,,™, + Di^,,m, (4.9) 
where we have introduced the following distributions with causal support: 

D = £)(.+)-£){-) £){-) 

m\,m2 mi m2 mi m2 

m\,m2 mi m2 mi m2 

^^^ri,m2 = ^S^^^'^'^^S - ^L^i^^'^'^^fc^- 

^^^.l,m2 = ^.^S^^^'^^S - ^^^^L^i^^^'^'^^fc^- (4.10) 

It is easy to prove that one can use the arbitrariness of the causal splitting of these distributions 
such that one can eliminate all anomalies of the Ward identities fl4.7|) - ()4.9|) . So it follows that 
only tree Feynman graphs can produce anomalies. 

We apply this strategy to the quantum gravity model from the preceding Section. 

First we consider the theory given by the interaction Lagrangian 

Theorem 4.1 The Lagrangian T = T^"^ gives a theory gauge invariant in the second order 
of the perturbation theory if we perform convenient finite renormalizations of the second-order 
chronological products 

T(x, y) = T^(x, y) + t Six - y)N{x) T^{x, y) = rj(x, y) + i 5{x - y)i N^{x) (4.11) 
where N and are some Wick polynomials. In particular: 

N = lQh'^''hP'^{dph,,){d,h) - 8/i^'^/i''^(9,V)(^"V) - 32/i'^^/i,,(9"/i'"^)(9^V) 

-lQh'"'h^p{d''h/){d^h) + lQuP{dpu'')u''{d,u^) 
+2m2 [hi^ - W^h^^K' + h^^^Kph':^h + W^hp^hP'^hp, - Ah^'h^php^h^''^ . (4.12) 

Proof: As it is known, the first step is to compute the causal commutator [T^{x),T{y)]. 
The anomalies are produced by two types of terms in T^{x) : (a) with the index /i appearing 
in a derivative 5^; (b) with the index /i appearing in the combination /i^p. Inspecting the 
expression from prop. 13.21 we have 

T, = T,{dph) + Tf{d,h^p) + mOpU,) + TlidpU,) + md,v,) 

+S''''{d,hpp) + SPh^p + --- (4.13) 

where by ■ ■ • we mean terms which do not produce anomalies and we have defined: 

Ti = -u^id^h) - {dPup)h + 2{dfu'')Kp 



T 



uxid^h"'^) + {d\xW^ - (9am")/i^^ - (9am^)/i^° 

= 2up{dPu'' 
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Sf = 4{daUx){d^h^^); (4.14) 

here we have imposed the symmetry condition 

rr^ip^a). (4.15) 



T2 = T2;«/9?7"^ (4.16) 



It is also convenient to denote 
with the explicit expression 

T2 = -2Ti. (4.17) 

One has to compute the commutator of d^h, d'^h"'^, d'^ii^, d^u^, d'^A^ with the 12 linear 
independent Wick monomials which appear in the expression of the total coupling T. Using 
the causal (anti)commutation relations from Section 4 we get : 

[d^h{x),T{y)] = i d^Dm{x - y)D,{y) + i d^d^D^ix - y)D'i{y) 
[d^H'^^ix), T{y)] = i d^D^ix - y)Df{y) + i d,d,D^{x - y)Df'^{y) 
[9>''(x), T{y)] = i d^Dmix - y)Dl{y) + i d^d,D^{x - y)D'^^''{y) 
[dX{x), T{y)\ = i d^D^ix - y)Dl{y) + i d^d,D^{x - y)D'^'''{y) 
[dX{x):T{y)] = i d^Dm{x - y)D',{y) + i d^d,Dm{x - y)D','''{y) (4.18) 



where 



where 



-A{d''u''){dpU^) + 4a {dfv"){dpv^) (4.19) 

Dl = -Ah^^idph) + Shpaid'h"") + 2^8" K) - ^Upid^u"). (4.20) 

Df = Vf + L>2?7°^, (4.21) 

Vf = -{d''h){d'^h) + 2{d"h''^){d^Kx) + 4.{dpN'''){d,hP^) 



+2{dPh''^){dph) - A{dPh''''){dphi) + 2(9X)('9V) + 2(9^Mp)(9" 



-4(a"wp) {df^vP) + 2m^ (^vP'^h'^p - (4.22) 



and 



D2 = -\{dph){dph) + {dPh''^){dpKp) - 2{d''h^^){dpKx) - 2{dpu''){dpu^) 



+2{dPv'') {dpv,) + (-h'^'^Kp + ^/i') , (4.23) 
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with 



D^3 = -4 



D'^'" = 4 



= Am{d''v'')ua 



One also has: 



with 



-2K'"{d,v")+mu''uP^. 



[h^p{x),T{y)] = i - y)Ep{y) + • 

= i df,d^Dm{x - y)Ep{y) + 



Ep = 2u''{daU 



p)- 



From these formulae we obtain 

= i 9^r>^(x - y)A{x, y) + i d^d,D^{x - y)A''{x, y) + ■ ■ ■ 

where 

5 

A{x,y) ^ T,{x)V,{y) + Tf{x)V2,aM + Y.^j{^)DrAy) + S''{^)E,{y) 

and 

5 

A\x,y) ^ T,{x)V'i{y) + T,,^p{x)Vf^\y) + Y.T,,p{x)Df\y) + 5^''(a;)E,(|/) 

j=3 
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where 

Vi = Di- 2D2 
Vl = Dl- 2D^. (4.38) 

From formula ()4.35|) we get the anomaly 

a{x, y) = 5{x - y)A{x) + [d,5{x - y)W{x, y) (4.39) 

By integration by parts we can rewrite this as follows: 

a(x, y) = 5{x - y)A{x) - dl [5{x - y)A^{y)\ (4.40) 

with 

A{x) = A{x) + dlA''{x,y)\y=, 

A''{x)=A^{x,x) (4.41) 

From the preceding formula we get 

5 

A = T,A + T^.e.pVf + Tr,pD^ + S^Ep + S^P{d,E,) (4.42) 

i=3 

where 

v^=v, + d,v'; j = i,2 

D,=Dj + d,D'^ J = 3, 4, 5, 6 (4.43) 

The total anomaly comes from the two commutators [T^{x),T{y)] + [T^{y),T{x)] and is 

A{x,y) = a{x,y) + a{y,x). (AAA) 

Now, the second-order gauge invariance condition (j2.45|) for n = 2 is fulfilled iff one can write 
this anomaly as follows: 

A{x, y) = 2zdQN{x, y) + d^N^^ix, y) + dlN>^{y, x) (4.45) 

where N{x, y) and N'^{x, y) are quasi- local Wick polynomials. One can show rather easily that 
this condition is equivalent to: 

A = idgN + df.N" (4.46) 

for some Wick polynomials (in one variable) and A^^ of ghost number and resp. 1; here A 
is given by the first formula (|4.4H) . From power counting considerations (j4.2j) one also has the 
limitations 

deg{N), degiN") < 6. (4.47) 
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The condition (j4.46j) is the basic condition which we will investigate from now on; it splits 
in distinct conditions for different sectors. 

a) We consider first the part Auhhh of the anomaly which is linear in and tri-linear in 
the field hap- From the expression of the anomaly we have 

Auhhh = + (5^M^)A + {d^u^)A,, (4.48) 

where 

A = -hVi + 2h''^Vf^p 



A^, = 2h^,V, - Ah-Pf)]^^^. (4.49) 



By the symbol hh as an index means that we consider only the bi-linear part in h^p of the 
corresponding expression. We want to find out if it is possible to write 

Auhhh = idQN^ + d^N^ (4.50) 

for some Wick polynomials A^^i and N'^ of ghost numbers and resp. 1. We will not consider 
the generic form of Ni and N'^ ; instead, we make the following ansatz 

iVf = f'" + (dpu,) t'^^'P (4.51) 

where the expressions t " are tri-linear in the fields h^p, have null ghost number and are limited 
by the power counting conditions 

degit^") < 5, deg{t^"'f) < 4; (4.52) 

we will also suppose that the expression t^^'^ does not contain terms with the factor rj^^f. 

If we consider in Ni only terms of the type hh{dh){dh) and hhhh, we have the following 
generic expressions of dgNi. 

I dgN, = {d'^un B^, + {d'd^un B^,, (4.53) 

where the expressions B... are tri-linear in the fields hap., have null ghost number and are limited 
by the power counting conditions 

deg{B^,) < 5, degiB^,,) < 4; (4.54) 

we can suppose that the expression B'^'^^ is symmetric in u, p. 

The explicit expressions for B,,, follow from the generic ansatz for A^i; they depend on some 
unknown coefficients which are to be determined. We substitute everything in the equation 
(j4.5Up and get the following system: 

A'^ = d.t"'' 
A"" + 7]"" A = B^"" + + dpt"'"' 

B^^P + Supt^'^P = (4.55) 



26 



From the last equation we can determine t'^'^'^ (if we suppose that it is symmetric in fi and 
p). We substitute in the second equation of the system and determine f^^. Finally we substitute 
this in the first equation of the system; it turns out that the system is consistent iff the following 
equation is true: 



(4.56) 



This is the basic equation which we now analyze. The expressions from the left hand side 
can be computed explicitly from the formulae ()4.49|) . The general strategy is to make an ansatz 
for A^i, compute the expressions and compute the right hand side of the equation. Then we 
get a system for the unknown coefficients of the finite renormalization A^^i. 

The ansatz for A^i can be guessed analyzing different contributions appearing in Auhhh and 
it is 

N, = /i/i'^'^/i'"^(9,V)(^-^) +/2^^'^""(^aV)(^"V) 

+fioh'"'KX-^^H + fuh'^'h^.hP^hp^ + fi2h'"'Kph^^hP''; (4.57) 

after performing the computation of dqNi we see that the expression Bfj,iyp does not have terms 
with the factor rj^p so tp^p = —Bypp does not have terms with the factor r/^^; this is consistent 
with the ansatz we have made for tp^p. The consistency equation (j4.56j) gives after long but 
straightforward computations: 



/i 



/2 



4 h = h 
' f 



-16 /5 = 16 /e 



f7 = -i 



-m 



10 



-m 
3 



11 



m 



'12 



-4m^. 



(4.58) 



b) We consider now the contribution Auuuh of the anomaly which is tri-linear in the ghost 
fields and linear in the field hajs- The explicit expression is: 



A 



uuuh 



{d^Up){d,3UP) + {daUp){dPUf3) + U^idadpUfs) 



X 



uxid'-u") \{dPd"K,)up + {d''K^){dpup) - (dphmd^u,) - h^^^idpO^u,) 



X 



-u^idxhP^) +ux{d^h''P) + {dxuP)h''^ 
We want to write this expression in the form 

Auuuh = i{dQN2) uuuh 

+ dpN^ 



(4.59) 



(4.60) 
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We make the ansatz 

N2 = fisU^idX^^'id^u,). (4.61) 
It is convenient to re- write the preceding formula (using "partial integration" ) as follows: 

Auuuh - i{dQN2)uuuh = d^N^ + KpU""" (4.62) 

where V^^ is tri-linear in the ghost fields. One can obtain after some computations that one 
can fix the constant /13 such that [/°^ = 0; namely we must have 

/i3 = 16. 

c) The contribution linear in the field of the anomaly is 



(4.63) 



Av = 8m 



(4.64) 



We have to subtract from this expression the contribution linear in from dqN^- The 
result is a divergence: 



d) Now we consider the contribution bi-linear in the field f^; we have 

= -8^(9^/1"^) + (9^Up)/i"^ - (^X)^""^ - {dy)hf"'] {d^v„){dpv'' 

One can write this expression as a total divergence: 

A 



(4.65) 



(4.66) 



(4.67) 



This finished the proof of the theorem. ■ 

We now consider the general case (i.e. 6 7^ 0). It is sufficient to consider the first two 
sectors. The sector uhhh modifies only the values of the coefficients fj. The anomaly Auhhh 
acquires the following expression: in (j4.48j) + (j4.49|) we have 



(4.68) 



where by ■ • ■ we mean the corresponding expressions from the case a = 1, 6 = multiplied by 
a^. If we perform the same computations as in the preceding theorem we obtain the new values 
of the coefficients fj: 

fi = Sa^ /2 = -4a2 f, = f, = -lQa' h = 16a^ f, = 8a^ = Sa' 

/s = (a^ + ^afe) h = - (o? + -afej w? fw = -(a^ + -abj rr? 



(4.69) 
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The sector uuuh gives the following expression of the anomaly: 

1 



where 



Auuuh = h rn^abAi + —m%^A2 

lb 



+2 



(4.70) 



(4.71) 



A2 = h 
-2/i^" 



u^{d\^)up - 2u^{dxUp)u^) - 2u^{dpUx)u^) 



(4.72) 



where by ■ • ■ we mean the expression (j4.59p multiplied by a^. If we try to write the total 
anomaly as a total divergence + coboundary, we obtain in the end: 6 = 0. 
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5 Conclusions 



The evidence for a dark energy in our Universe motivates the introduction of the cosmological 
constant A into Einstein's equations. The corresponding Einstein - Hilbert Lagrangian is given 

by 



- — v^(i?-2A), K2 = 327rG 



(5.1) 



where G is Newton's constant and g = det{g^i^). We want to expand this in powers of n. For 
this purpose it is convenient to use the so-called Goldberg variables 



Now we write this metric tensor as 



(5.2) 
(5.3) 



(5.4) 



Then we can expand the determinant g and the Ricci scalar R in (1) in powers of details 
can be found in (SHI, Sect. 5. 5. 

The first four orders in i.e. 0(/€~^, . . . 0(/€^) come out to be: 

Lb = ^ + -{d^h - dadf^h'^" + Ah) 



K 



+4A 



t/37' 



L32V"^ "^"'^y ■ 12' " 32 "'^ ' 4!16 8 

Here the terms L^^^ and L*^^) without A are the same as in the ordinary massless gravity: 

L(^) = -hi'^l^{d^h){dph) + ^/i'^^(9^/i"^)(9^/i"^) + /i"'^(9^/i°^)((9^/i^^) 



(5.5) 



(5.6) 



L(2) = -h^,h':f,{d,h^nidf,h''n - -hpph''^{d^h'^){d''h) 



(5.7) 
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The terms 0(/t ^) and 0(k ^) in the Lagrangian (|5.5|) are irrelevant because they give no 
equation of motion. The term 0{kP) and quadratic in h 

^ l(d^'h-P){d,Kf,) - {d,h"^){df,h^J - ^{d^h){d''h) + A(y - h^^h^f,) (5.8) 
gives the following Euler-Lagrange equation 

d^h"'^ - d^{d^h^^ + d'^h^^') - ^r]'"^d^h - A(r7"^/i - 2/i°^) = 0. (5.9) 
Taking the trace we find 

d^d^h""^ = -\d^h - Ah. (5.10) 

Differentiating ()5.9|) by da and substituting ()5.10|) we derive the Hilbert gauge condition 

dah'''^ = 0. (5.11) 

Then ()5.1Up reduces to the Klein-Gordon equation 

d'^h + 2Ah = (5.12) 

and from ()5.9p we obtain the Klein-Gordon equation for the tensor field 

a^/i"^ + 2A/i"^ = (5.13) 

This means the graviton becomes massive with mass 

= 2A. (5.14) 

Taking the current value of A one finds a tiny mass (~ 10~^^ eV) for the graviton, so that 
our massive theory passes all direct tests of general relativity. However, gravitational radiation 
requires a thorough investigation. 

The cubic part O(k^) in ()5.5p gives the coupling. Using ()5.14p this agrees exactly with the 
pure graviton coupling terms in ()3.61|) . if we multiply with an overall factor —4. Similarly, 
the pure graviton couplings in the quartic part O(k^) agree with ()4.12p without the ghost 
terms, if we multiply by 16. This shows that our massive gravity is the quantum gauge theory 
corresponding to classical gravity with a cosmological term. Of course, the ghost couplings 
cannot be derived from the simple Lagrangian above. To guess the right Lagrangian with ghost 
fields (fermionic and bosonic) is not quite straightforward. The best way to deduce the full 
theory is quantum gauge invariance, instead of some classical heuristics. 

In spin 1 gauge theories with massive gauge fields a Higgs field is necessary to "generate" 
the masses. In the framework of quantum gauge theory it is simply needed to restore gauge 
invariance to second order [21]; first order gauge invariance holds without the Higgs couplings. 
It was a surprise for us when we found that massive gravity is gauge invariant to first and 
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second order without a Higgs field. In a way this is even disappointing because a gravitational 
Higgs field would be a nice candidate for the non-barionic dark matter. 

A last remark is concerned with the mass zero limit of our theory. As noticed above, in 
the limit m ^ the bosonic ghost does not completely decouple from the graviton. In 
fact, the term 4:h^'^{d^Vp){dyV'') survives in ()3.61|) . That means the resulting massless theory is 
not identical with usual quantum gravity as discussed in [21], for example. This leads to the 
conclusion that there exists at least two different quantum gauge theories which correspond to 
classical (massless) general relativity, one with an additional Bose field and one without. 
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